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Abstract 

New research is presented, and previous research is reviewed, on the emission and 
absorption of interstellar aromatic hydrocarbons. Emission from aromatic hydrocarbons 
dominate the mid-infrared emission of many galaxies, including our own Milky Way galaxy. 
Only recently have aromatic hydrocarbons been observed in absorption in the interstellar 
medium, along lines of sight with high column densities of interstellar gas and dust. Much 
work on interstellar aromatics has been done, with astronomical observations and laboratory 
and theoretical astrochemistry. In many cases the predictions of laboratory and theoretical 
work are confirmed by astronomical observations, but in other cases clear discrepancies 
exist which provide problems to be solved by a combination of astronomical observations, 
laboratory studies, and theoretical studies. The emphasis of this paper will be on current 
outstanding puzzles concerning aromatic hydrocarbons which require further laboratory 
and theoretical astrochemistry to resolve. This paper will also touch on related topics 
where laboratory and theoretical astrochemistry studies are needed to explain astrophysical 
observations, such as a possible absorption feature due to interstellar "diamonds" and the 
search for fullerenes in space. 
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1. Introduction 

Interstellar atomic gas, molecular gas, and dust grains, where grains are solid particles 
with sizes of 1 - 1000 nm, pervade the space in between the stars known as the interstellar 
medium (ISM) . The grains and molecules in the lower density regions of our Galaxy, called 
the diffuse ISM, are believed to form primarily in the outflowing gas from dying stars with 
masses similar to that of our own Sun. The final stage of this type of stellar death is a shell 
of ejected gas, known as a planetary nebula, which surrounds the stellar corpse, known as a 
white dwarf. Grains and molecules find it difficult to form in the general ISM, which has 
typical gas densities of 1 - 10 atoms cm~^, and thus extremely low chemical reaction rates. 
The outflows from dying stars have much higher densities (> 10^ atoms cm^^), which make 
chemical reactions leading to formation of molecules and grains much more probable. 

In our universe, the pattern of cosmic elemental abundances is determined first by 
the hydrogen and helium formed in the Big Bang, when our universe began, and secondly 
by nuclear reactions within stars, which account for virtually all the other elements. As a 
result of the pattern of common nuclear reactions within stars, only a very restricted set 
of elements are important in interstellar chemistry. Figure 1 shows the cosmic abundance, 
relative to hydrogen, of all elements which have an abundance larger than 10"^ relative to 
hydrogen [1]. It can be seen that the most important elements for interstellar chemistry are 
H, O, G, N, Mg, Si, Fe, S, Al, Ga, Na, and Ni. The inert gases He, Ne, and Ar are also 
relatively abundant but do not interact chemically. Other elements exist in the universe, 
but their scarcity leads to their playing a minor role in interstellar chemistry. By far the 
most abundant interactive elements are H, O, G, and N, and thus it is no surprise that the 
overwhelming majority of interstellar molecules detected to date are primarily composed of 
these elements. 

The interstellar infrared emission features (lEFs), at 3.3, 6.2, 7.7, 8.6, 11.3, and 12.7 
lim (790, 890, 1160, 1300, 1610, and 3040 cm~^), dominate the mid-infrared emission of our 
own and other galaxies. Figure 2 shows the 600 - 2000 cm~^ spectrum of one region within 
our own Galaxy, NGG 7023 [2] . This figure shows most of the lEFs. Note that this figure 
is plotted versus cm~^, to be helpful to chemists, but that the astronomical convention is 
to plot spectra versus /im and that the astronomical convention will be followed for the 
remainder of this paper. 

The lEFs, which have widths broader than natural broadening mechanisms for fines 
in atomic gas or simple molecular gas, were mainly first discovered by astronomical 
observations in the 1970's [3-6], with the 12.7 /xm IFF discovered only after higher 
spectral resolution observations [7] were able to distinguish it from Ne"*" emission at 12.8 
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lim (780 cm~^). The lEFs have been referred to by a variety of names in the hterature, 
with "feature" and "band" used interchangeably: the "Tinidcntificd emission features," 
the "unidentified infrared bands," the "infrared emission features," the "aromatic infrared 
bands," and the "polycychc aromatic hydrocarbon (PAH) bands." The origin of the last 
two names will become apparent in later paragraphs. 

The lEFs are believed to arise in either interstellar grains or interstellar molecules, 
and because their emission can account for up to 30 - 40% of the overall emission of 
galaxies, solving problems related to their composition, size, and origin are essential for 
understanding the overall energy budget of radiation from galaxies. Within our galaxy, the 
lEFs are observed in a wide variety of environments. One important clue to the composition 
of the carrier of the lEFs, however, is the astronomical observation that the strength of 
the 7.7 /im lEF shows a strong correlation with the ratio of carbon to oxygen in planetary 
nebulae [8 - 10]. Ordinarily there is more oxygen than carbon in interstellar gas and in 
stars, which results in an interstellar chemistry in which common elements most often 
combine with hydrogen and oxygen to form molecules and grains. The heavier elements 
in such an oxygen-rich environment generally condense into solid dust grains composed 
of oxides of these elements, such as silicate grains which are commonly observed in the 
ISM. In some dying stars, however, convection within the star dredges up so much carbon, 
produced in nuclear reactions at the star's center, that the composition of the outflowing 
gas from a dying star changes to a peculiar situation in which the gas has more carbon 
than oxygen. This carbon-rich gas alters the chemistry of grain and molecule formation in 
the outflowing gas, preferentially forming molecules and grains which are rich in hydrogen 
and carbon. The observed correlation between the 7.7 /xm lEF strength and the C/0 ratio 
in planetary nebulae strongly suggests that the lEFs form from dust and/or molecules 
condensed in the outflows from dying stars with carbon-rich outflows. Once the lEFs, like 
other grains and molecules, are formed, the gas ejected from the dying star slowly mixes 
with the surrounding gas and dust within the galaxy and becomes part of the general 
ISM. The ISM in galaxies, therefore, contains a mixture of molecules and grains formed in 
oxygen-rich environments and carbon-rich environments. 

Astronomical observations show strong emission from the lEFs in many regions of 
star formation, where dust and gas are interacting with the ultraviolet (UV) and visible 
radiation emanating from young stars or from stars in the process of forming (protostars) . 
Two examples of star formation regions with strong lEF emission are reflection nebulae, 
where grains within a cloud of gas and dust scatters the light of a nearby young star, and 
photodissociation regions surrounding H II regions (see reviews by [11, 12]). An H II region 
is the astronomical name for a cloud of ionized hydrogen gas surrounding a young star hot 
enough to ionize hydrogen; a photodissociation region is the astronomical name for the 
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surface of a cloud of molecular hydrogen gas, where the UV radiation from a hot star is in 
the process of photodissociating the molecular hydrogen. Regions of the ISM are commonly 
referred to by the state of their hydrogen gas - neutral (H I), ionized (H II), or molecular 
(H2) - because 71% of the mass of interstellar gas is composed of hydrogen. The hydrogen 
gas in reflection nebulae can be either neutral or molecular, but not ionized. 

Astronomical mapping and images of the lEFs in the H II regions and adjacent 
photodissociation regions in the star-forming Orion Nebula and the planetary nebula NGC 
7027 suggest that the lEFs are concentrated in a narrow zone between the ionization 
front and the photodissociation front [13, 14]. The ionization front is where ionized 
hydrogen abruptly changes to neutral hydrogen at a distance from the ionizing star where 
all of the ionizing photons from the star has been absorbed by intervening hydrogen gas 
and re-radiated as several lower energy, non-ionizing photons. The ionization front can 
be delineated by astronomical observations of recombination lines of neutral hydrogen; 
such lines are strong in ionized hydrogen but absent from neutral hydrogen. The 
photodissociation front can be traced by astronomical observations of II2 emission; the 
same UV photons which can dissociate H2 have an alternate path in which they cascade 
through the electronic levels of H2 to populate vibrational-rotational levels of H2, which are 
then observed as H2 emission lines in the near-infrared (1 - 2.5 /xm or 4000 - 10,000 cm~^). 

The lEFs are also observed in regions of low UV intensity, such as the surfaces of 
molecular clouds and the general ISM where the illumination source is the diffuse interstellar 
radiation field [15 - 17]. The ubiquity of lEF emission within our Galaxy, within star 
formation regions, within carbon-rich planetary nebulae, and within the general ISM, 
makes understanding their carriers an essential part of understanding how our Galaxy 
emits, absorbs, and reprocesses radiation. The prevalence of lEF emission in star formation 
regions, and the recent detection of lEF absorption toward protostars [18 - 23], also 
illuminates the important role that the lEF carriers play in the material from which stars 
and planets form. 

The lEFs were first proposed to be due to aromatic hydrocarbons by Duley & Williams 
[24] , based on wavelength coincidences between the lEFs and laboratory spectra of aromatic 
hydrocarbons. Another clue to the composition of the lEF emitters was contributed by 
Sellgren et al. [25] and Sellgren [26], when we discovered 1-5 fim. (2000 - 10,000 cm~^) 
continuum emission in refiection nebula that was not due to scattered starlight. Starlight 
which is scattered from grains is strongly polarized, and this infrared continuum emission in 
refiection nebulae was unpolarized [27]. This continuum emission had a color temperature 
around 1000 K, and was spatially associated with lEF emission in these sources. The color 
temperature of the continuum emission showed no decrease with distance from the exciting 
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star, as would be expected for thermal emission by dust in equilibrium with the star's 
radiation field. Furthermore the stellar radiation field was far too weak to heat grains to 
equilibrium temperatures of 1000 K at the observed nebular distances. Temperatures of ~ 
30 - 60 K were expected instead, and in fact dust in this temperature range had already 
been observed [28, 29] at far infrared wavelengths (40 - 400 fira or 25 - 250 cm"-*^). We 
proposed a new paradigm for interstellar dust, in which the dust size distribution contained 
a component of 1 nm sized grains which were so small that a single UV photon absorbed 
from a nearby star could briefiy heat the grain to temperatures as high as 1000 K. Tiny 
particles of this size has been proposed earlier [30, 31]. These tiny grains spend the vast 
majority of their time at low temperature, but when they undergo stochastic heating they 
can for a microsecond or millisecond radiate emission at high temperatures [25, 26]. This 
model explains the high color temperature of the continuum emission, the constancy of the 
color temperature with distance from the exciting star, and the lack of polarization. Once 
we had proposed this model, several groups quickly realized the importance of this emission 
mechanism to the newly emerging data from the Infrared Astronomical Satellite (IRAS), 
launched in 1984, because the tiny grains, as they cooled radiatively, would be characterized 
by progressively lower color temperatures with time and thus contribute significantly to 
emission in the IRAS 12 and 25 //m (400 and 830 cm~^) broad-band photometric filters [32 
-34]. 

The real breakthrough in understanding the lEFs, however, was made by Leger & 
Puget [35] and AUamandola et al. [36]. They were the first to combine the aromatic 
hydrocarbon idea [24] with a 1 nm particle size [26] , and to propose that the lEFs were due 
to polycyclic aromatic hydrocarbon (PAH) molecules with a size of roughly 1 nm or smaller. 
Other aromatic materials have been proposed for the lEFs, including tiny grains composed 
of hydrogenated amorphous carbon, quenched carbonaceous composite, and coal-like 
materials [37 - 42]. Virtually all astronomers working on the lEFs agree on the aromatic 
nature of the lEFs. There is still debate over the exact nature of the interstellar aromatics, 
but many favor the PAH hypothesis because it so neatly ties together well-studied aromatic 
materials with the 1 nm sizes suggested by the continuum emission associated with the 
lEFs. 

The goal of the astronomical observational research of myself and my collaborators 
has been to test specific predictions of the PAH hypothesis, and to urge chemical theorists 
and laboratory astrochemists to make modifications in the PAH hypothesis when the 
predictions and observations do not match. This interplay between the observations, theory, 
and laboratory work has proven very fruitful in refining our understanding of the size, 
composition, ionization state, and hydrogenation of the material responsible for the lEFs in 
the ISM. 
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1.1. A Challenge: Approximating Interstellcir Conditions in the Laboratory 



One of the most difficult challenges in matching astronomical observations and 
laboratory measurements of candidate materials for the lEFs is the vast gulf between 
interstellar conditions and laboratory conditions. The densities in interstellar space, even 
in the densest interstellar clouds, are so low that they exceed the capabilities of the best 
vacuum systems on Earth. The vast majority of published laboratory work has been 
measurements of solid-phase, room temperature aromatic materials in absorption, which 
arc then compared to astrophysical observations of ~1000 K aromatic materials in emission 
in extremely low density environments. Some progress has been made in approximating 
the low densities of interstellar conditions by use of inert-gas matrix isolation techniques, 
but even in these cases matrix interactions introduce frequency shifts in the wavelengths of 
absorption lines of unknown magnitude and direction. Other progress has been made in 
measuring aromatic molecules in emission in the laboratory, which has helped to quantify 
the temperature dependence of wavelength shifts in PAHs. Recent evidence, discussed 
below, suggests that the PAHs that have been best-studied in the laboratory are too small 
to approximate interstellar PAHs. Furthermore, astrophysical theoretical calculations 
predict that PAHs should be ionized in the best-observed regions of the ISM. Both chemical 
quantum calculations and laboratory measurements demonstrate marked changes in the 
relative intensities and central wavelengths of PAH bands with ionization state. All of these 
mismatches between interstellar conditions and laboratory conditions make it difficult to 
achieve a firm identification of the material(s) responsible for the lEFs. 



2. The lEFs and the Hardness of the UV Radiation Field 



One of the early predictions of the PAH model is that PAH molecules should only be 
excited by UV radiation, since laboratory absorption curves of small, neutral PAHs show 
a sharp cutoff in their absorption at UV wavelengths with little or no absorption at visible 
wavelengths. Scllgren et al. [43] tested this prediction by using observations from the IRAS 
satellite to search for infrared emission in reflection nebulae excited by stars of different 
effective temperature. Teg. We first verified that each reflection nebula was indeed excited 
by the central star of the nebula, by requiring that the dust temperature derived from the 
the ratio of 60 iim (170 cm~^) intensity to 100 iim (100 cm~^) intensity reach a peak value 
at the star. The emission at 60 and 100 /im is thought to primarily come from grains, 
with a size around 100 nm, which are in equilibrium with the radiation field, and whose 
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temperature therefore will increase with proximity to the heating star. We then examined 
the 12 //m emission of each reflection nebula, at a nebular position offset from the star 
so that stellar emission would not contaminate the 12 /xm data, and compared it to the 
total infrared emission at the same spatial location. For each nebula we then measured 
i?(12/total), the ratio of the total flux in the IRAS 12 ^m band to the total infrared flux. 
The IRAS 12 yum band we assumed to be dominated by lEF emission and its associated 
continuum emission, as had been shown by mid-infrared spectroscopy for some of the 
sources in our sample [34]. We were expecting i?(12/total) to show a precipitous drop for 
cool stars {T^s < 10,000 K), where the fraction of total stellar radiation radiated in the UV 
was small. Much to our surprise, the value of i?(12/total) was independent of Tes, over 
the range T^s = 5,000 K - 22,000 K, as illustrated in Figure 3. Our observational results 
required that the material responsible for emitting the IRAS 12 yum emission in reflection 
nebulae absorb not only at UV but also at visible wavelengths. This discovery altered the 
astronomical community's perception of PAHs as being small in size and neutrally charged, 
and drove theorists and laboratory astrochemists to consider both larger PAHs and ionized 
PAHs in their models and laboratory work. This is because increasing the PAH size and 
ionizing PAHs both have the effect of extending the absorption cross-section of PAHs out 
to visible wavelengths. 

The 12 iim IRAS filter bandpass is wide and encompasses several lEFs as well as 
their associated continuum. The material responsible for the continuum has never been 
identified, and the possibility exists that it is not due to the same material that produces 
the lEFs. In this case, it might be possible that the IRAS 12 //m emission we observed in 
refiection nebulae illuminated by cool stars is due to continuum emission alone without 
any contribution from the lEFs. Thus, when the Infrared Space Observatory (ISO) was 
launched in 1995, we embarked on an observational program to obtain spectra of some 
of the reflection nebulae studied by Sellgren et al. [43], particularly reflection nebulae 
illuminated by both cool and hot stars which had similar values of i?(12/total) . These 
observations have been made with ISO's mid- infrared camera (ISOCAM) combined with 
its circular variable filter (CVF), which allow us to obtain low-resolution (A/ A A = 40) 
spectra, at 5 - 15 //m (600 - 2000 cm~^), simultaneously at roughly a thousand spatial 
positions across each nebula. The ISO spectra of these refiection nebulae, therefore, should 
unambiguously determine whether the lEFs are present in these sources. 

The first ISO results from Uchida et al. [44] was our discovery that lEFs are clearly 
detected in vdB 133, a refiection nebula illuminated by a binary system with very little 
UV radiation. The binary system, a luminous star with T^s = 6,800 K plus a fainter star 
with Tcff = 12,000 K, provides a ratio of UV (A < 400 nm) to total stellar flux which is a 
factor of four lower than more typical refiection nebulae which are illuminated by stars with 
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Teff ~ 20,000 K. Yet, despite the softer radiation environment, the lEF spectrum in this 
UV-poor environment is very similar to lEF spectra observed in sources with much harsher 
UV environments, as shown in Figure 4. 



2.1. New Results: Laboratory Analogs and the IRAS Data on R(12/total) 



I am currently collaborating with J. Pizagno, K. Uchida, and M. Werner on a 
comparison of various laboratory and theoretical candidates for the lEF carriers with the 
IRAS observations [43] of i?(12/total), to quantify which materials can reproduce the 
lack of dependence of i?(12/total) on T^s (Pizagno et al. 2000, in preparation). We are 
convolving the UV and visible laboratory and theoretical absorption curves of different 
materials with the energy distributions emitted by stars of different Teff, and comparing the 
results to our IRAS observations. Figure 3 shows one sample result of our calculations, 
where we compare the astronomical observations to the predicted results for two small PAH 
molecules, neutral naphthalene [45] and singly ionized naphthalene (from F. Salama 1999, 
private communication). These new results show that small neutral PAHs are completely 
unable to provide enough visible absorption to explain the lEF emission observed around 
UV-poor sources. They also show that while ionized PAHs have more visible absorption 
than their neutral counterparts, small PAHs even when ionized also cannot explain the 
IRAS observations. 



2.2. Need for Data and Theory for Larger PAHs and Varying Ionization 



It seems likely that a combination of ionization and larger PAH size will be required to 
explain the astronomical observations. We strongly encourage laboratory astrochemists to 
measure parameters for larger PAHs, both neutral and ionized. We also strongly encourage 
theoretical astrochemists to calculate parameters for larger PAHs, both neutral and ionized, 
particularly for PAH sizes and/or ionization states that are not currently accessible by 
laboratory techniques. 



3. PAH Ionization 
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Uchida et al. [46] have just completed a more comprehensive ISO study of lEFs in 
reflection nebulae, in which we have made Tinexpected observations which have important 
ramifications for models of PAH ionization. PAH ionization models predict that PAHs 
should be primarily positively charged in regions of high UV radiation, such as refiection 
nebulae, and primarily neutral or negatively charged in regions of low UV radiation, such 
as the diffuse ISM and cirrus clouds [47 - 50]. The UV radiation intensity is generally 
characterized by the variable Gq, where Go = 1 corresponds to the UV radiation field in 
the solar neighborhood. Laboratory and theoretical studies show that the ratio of PAH 
emission at 6 - 10 //m (1000 - 1700 cm"^) to the PAH emission at 10 - 14 //m (700 - 1000 
cm~^) is a sensitive function of PAH ionization state, with this ratio varying by at least a 
factor of 10 between ionized and neutral PAHs [51 - 59]. 

We have measured the ratio of lEF emission at 6 - 10 iim to the lEF emission at 
10-14 //m, observed in refiection nebulae with Go varying from 20 to 6x10^. Figure 5 
shows that we see little or no evidence for any change in this ratio with Go (and thus PAH 
ionization state). This provides a major puzzle for the PAH hypothesis. 



3.1. PAH Ionization: Astrochemical Theory and Laboratory Data Needed 

There are several ways to resolve the discrepancy between the observations and theory 
for the predicted dependence of the PAH ionization state on the UV field. The calculated 
ionization states of PAHs are somewhat sensitive to the assumed PAH size. The amount 
by which the ratio of fiuxes at 6 - 10 //m to 10 - 14 //m changes also depends on PAH 
size, with smaller changes in the ratio for larger PAHs. The calculated ionization states 
of PAHs are also very sensitive to the uncertain PAH recombination rate, which requires 
better laboratory measurements and/or better theoretical predictions for the recombination 
rates, particularly as a function of PAH size. These are all areas in which astrochemistry in 
the laboratory and theoretical calculations can really make a difference in interpreting the 
astronomical observations. 



4. The Full- Width at Half Maximum of the 7.7 //m lEF 

Another unexpected discovery Uchida et al. [46] have made with our recent ISO 
observations concerns spectral changes in the lEFs at low levels of illumination. We have 
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quantitatively compared the lEF spectra (Fig. 4) of three sources we observed ourselves, 
vdB 133, vdB 17 (NGC 1333), and vdB 59 (NGC 2068), to the lEF spectra of pubhshed 
spectra of NGC 7023 ([2]; Fig. 1) and p Oph [15]. Wc have found no evidence for any 
systematic spectral differences with Teg. Our new ISO observations, however, find that 
the full width at half maximum (FWHM) of the 7.7 pirn. lEF is dependent on the distance 
between star and nebula in vdB 17 [46]. Figures 6 and 7 show that in the most distant 
regions of vdB 17, corresponding to the lowest UV illumination levels, Gq = 20 - 60, the 
7.7 //m lEF becomes significantly broader and begins to blend with the 8.6 //m lEF. This 
effect has now been observed in a second reflection nebula, Ced 201 [60], although they 
do not give the Go values at which the broadening of the 7.7 //m lEF becomes significant. 
None of the other lEFs appear to change width with stellar distance in either source. We 
also find no broadening of the 7.7 pirn. lEF or other lEFs in vdB 59, down to Gq = 200. 
More significantly, the spectrum of p Oph [15] at Go = 40 shows no sign of a broader 7.7 
Hm lEF or any blending of the 7.7 /xm lEF with the 8.6 fim. 

The unexpected broadening of the 7.7 /xm lEF at low UV illumination levels, in some 
sources but not others, is very much a mystery. One possibility is that the intermittent 
broadening of the 7.7 p,m feature depends on two parameters, one being Go, for instance, 
and the other being some factor such as Tes- Stellar T^s seems like a good starting point 
for investigation, because vdB 17 and Ced 201, which are illuminated by stars with T^s = 
10,000 - 11,000 K, show the broadening of the 7.7 //m lEF at low Gq values, while p Oph, 
which shows no such broadening at Gq = 40, is illuminated by a pair of stars with T^g — 
22,000 K. 

Two other ideas seem like promising avenues to explore. One is the recent ISO 
observations of Moutou et al. [61], illustrated in Figure 8. We find that the 7.7 iim lEF, 
when observed at higher spectral resolution (A/ A A = 1800), appears to break up into at 
least three different features [61], at 7.45, 7.6, and 7.8 /xm (1280, 1320, and 1340 cm-^). 

Previous astronomical spectra at A/AA = 100 - 200 [62, 63] resolved the 7.7 //m lEF into 
two features at 7.6 and 7.8 p,m. Other astronomical observations [9, 10] at a lower resolution 
of A/AA = 50, where only the central wavelength of the 7.7 p,m lEF could be measured, 
demonstrated that the 7.7 p.m lEF central wavelength depends on the physical conditions 
in the emitting region (whether or not the hydrogen gas is ionized, for instance). Since the 
results of [9, 10] almost surely result from subtle changes in the relative strengths of the 
subcomponents of the 7.7 //m lEF, then the broadening of the 7.7 iim lEF observed at low 
resolution [46] could reflect more obvious changes in the relative strengths of the different 
subcomponents observed at higher resolution [61 - 63]. 

The other idea deserving further exploration is related to the fundamental reason why 
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the lEFs are so much broader than gas molecules or grains moving at the Doppler speeds (1 
km or less) typical of these interstellar regions. There is no clear consensus on whether 
the lEF widths are due to physical processes (such as timescales for internal molecular 
processes, which depend on the size and temperature of the molecule) or to mixtures of 
aromatics with slightly different compositions, sizes, and central wavelengths. Observations 
of the lEFs in absorption suggest that the mixture of lEF carriers may be dominant in 
determining the central wavelength and width of the lEF emission, as discussed in the next 
section, and if this is the case, then the astronomical observations of changes in the central 
wavelength or FWHM of the 7.7 //m lEF could reflect changes in the size distribution or 
composition of interstellar aromatics under different physical conditions. 



4.1. The 7.7 Micron lEF Width: Need for Laboratory and Theoretical Work 

Astrochemical work is desperately needed to formulate some theoretical framework 
for why the 7.7 fim width changes, so that astronomical observations can be made to test 
this framework. As described above, the intermittent broadening of the 7.7 fim lEF could 
depend on two parameters, one being Go, for instance, and the other being some factor 
such as stellar temperature. Can laboratory or theoretical work provide any insights into 
why these two parameters, or some other set of parameters, might affect the 7.7 /im width? 
Does the broadening of the 7.7 iim lEF observed at low resolution reflect changes in the 
relative strengths of the different subcomponents at 7.45, 7.6, and 7.8 //m observed at 
higher resolution [61]? Or do the relative strengths of the different subcomponents of the 
7.7 fim lEF stay the same, but the width of each subcomponent change? What physical 
change would be needed in the lEF carriers for the width to change? Could it be that 
the size distribution or composition of interstellar aromatics change with different physical 
conditions? Or is it some other phenomenon entirely? The answer to the questions posed 
by observations of the 7.7 //m lEF are rooted in the ongoing debate as to the origin of the 
lEF widths themselves. 



5. Aromatics in Absorption 

Until recently, all astronomical observations of aromatic hydrocarbons in the ISM were 
of the lEFs in emission. The absorption signatures of the lEFs had not been seen. We 
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have begun a program to try to trace the abundance and chemical changes of aromatic 
hydrocarbons during their residence time in molecular clouds. A broad 3.25 fim (3080 
cm^^) absorption feature along the line-of-sight to the protostar Mon R2/IRS-3 was 
tentatively detected, and its discovery confirmed with better spectra, by Sellgren et al. [18, 
19]. We proposed that this feature is due to the C-H stretch absorption of cold aromatic 
hydrocarbons within the molecular cloud. Laboratory work on PAHs [64] has shown that 
the central wavelength of PAH features shifts to increasing wavelength with increasing 
temperature. This predicts that aromatics in absorption should be blue-shifted by about 
the observed amount, relative to aromatics in emission, because aromatics in emission are 
emitting at temperatures of ~1000 K, while aromatics in absorption are absorbing at much 
lower temperatures, '^10-80 K. In Brooke et al. [21, 22] we have now detected the 3.25 
/im absorption feature toward a total of five protostars embedded in molecular clouds. 



5.1. New Results: Compciring Aromatics in Emission and Absorption 

The current observations of aromatics in absorption, compared to aromatics in 
emission, at first glance seem to have a simple interpretation. But there is an interesting 
complication. Ground-based [65] and ISO [66] spectra of sources near the Galactic Center 
reveal a narrow absorption feature near 3.28 /im (3050 cm~^), which these authors identify 
with aromatic hydrocarbons. If this feature is due to the 3.29 iim (3040 cm~^) lEF emitter, 
observed in absorption through the diffuse ISM in our Galaxy, it should spend the vast 
majority of its time at cold temperatures and thus have a wavelength blue-shifted relative 
to the 3.29 fim lEF in emission, similar to the blue-shifted wavelength (3.25 /im; 3080 
cm~^) observed in absorption toward molecular clouds. The longer wavelength of the 3.28 
lim absorption suggests either that the aromatics are at a higher temperature, perhaps 
due to localized heating by sources in the Galactic Center, or that the composition, sizes, 
or ionization state of the aromatics observed in the diffuse ISM are different from the 
aromatics observed in molecular clouds. 

The differences in observed width are also confusing. Laboratory studies [64] have 
shown that as individual PAHs are heated, their central wavelengths both shift to longer 
wavelengths and the width of the PAH band broadens. Figure 9, which has not been 
published elsewhere, illustrates the contrast between the short wavelength and broad width 
of the 3.25 fj,m aromatic absorption in molecular clouds [19], and the longer wavelength and 
narrow width of the 3.28 /im aromatic absorption towards the Galactic Center [66]. Both 
spectra have A/AA 1000. Figure 9 also shows a spectrum of the 3.29 /xm lEF in emission 
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in the reflection nebula NGC 7023, from the work in progress of myself, C. Moutou, A. 
Leger, L. Verstraete, M. Werner, M. Giard, and D. Rouan. This spcctiTim, with A/AA = 
1000, demonstrates that the observed width of the 3.29 /im lEF in emission (Sellgren et 
al. 2000, in preparation) is similar to or narrower than the 3.25 fim feature in absorption 
[19], again in contradiction to the temperature model. This suggests that the observed 
wavelength shifts and widths are not due to temperature efi'ects but rather to composition, 
size distribution, ionization, or other effects. 

This mystery is underscored by recent observations of an absorption feature at 6.18 
fim (1620 cm~^), believed to be the aromatic C-C stretch, toward five hot, mass-losing 
stars and two Galactic Center sources in ISO spectra [67]. A possible absorption feature 
at 6.24 //m (1600 cm~^) toward several protostars has also been tentatively identified as 
the C-C stretching mode of aromatic hydrocarbons [20, 23). If the 6.18 iim absorption 
toward the sources observed by Schutte et al. [67] arises from aromatics in the diffuse ISM, 
the aromatics should again be cold with blue-shifted, narrow absorption features. The 
6.24 //m absorption feature observed by Keane et al. [23] toward protostars in molecular 
clouds, if due to aromatics, should also arise from cold aromatics with blue-shifted, narrow 
absorption features. The astronomical observations, however, do not fit this picture. The 
6.18 /xm absorption feature observed through the diffuse ISM [67] is blue-shifted relative to 
the 6.22 /im (1610 cm~^) lEF in emission [68], as expected, but the 6.24 /xm absorption in 
molecular clouds is observed at a similar or shghtly longer wavelength as the 6.22 //m lEF 
in emission. Furthermore, the 6.18 /im absorption feature observed through the diffuse ISM 
is significantly narrower than the 6.24 /im absorption feature observed towards molecular 
clouds, as shown in Figure 10 of Keane et al. [23]. Again this is strong evidence that both 
the central wavelengths and widths of aromatics observed in absorption and emission in 
the ISM arc determined not by temperature or molecular physics of a single molecule, but 
rather by the composition mixture, size distribution, or mix of ionization states of aromatics 
observed along a particular line-of-sight. 



5.2. Need for Laboratory Work on PAHs in Water Ice Matrices 

One possible way to reconcile the temperature picture of the wavelength and width 
of aromatic absorption and emission with the astronomical observations is to consider 
whether aromatics in molecular clouds might be embedded in ice mantles on grains. Solid 
H2O ice is observed toward all sources in molecular clouds with detected 3.3 or 6.2 //m 
aromatic absorption features. New astrochemical laboratory measurements at NASA's 
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Ames Research Center are currently underway to measure the effect of a surrounding ice 
matrix on PAH absorption (M. Bernstein 2000, private comnnmication). These results, 
once published, may clarify our understanding of aromatic absorption in molecular clouds. 
Measurements are needed for both neutral and ionized PAHs. 



6. Interstellcir "Diamond- like" Ccirbon 

The broad 3.47 /xm (2880 cm~^) absorption feature (FWHM pa 0.1 /xm or 80 cm~^) 
was first noted in ground-based spectra of four protostars by AUamandola et al. [69] . They 
suggested that the feature might be due to the C-H stretch absorption of solo hydrogens 
attached to sp^ bonded carbon clusters, the "diamond" -hke form of carbon. The feature 
was present in every molecular cloud source looked at by Brooke et al. [21]. Interstellar 
microdiamonds, formed outside our solar system, have been identified in meteorites [70, 71] 
and so must therefore exist and survive in the diffuse ISM. Our observations [21, 22] of the 
3.47 ^m. absorption feature, however, reached the startling conclusion that the optical depth 
of the 3.47 /xm absorption feature is not correlated with the depth of silicate absorption, as 
might be expected for two refractory minerals such as diamonds and silicates (Figure 10). 
Instead, we found that the 3.47 /xm optical depth is strongly correlated with the depth of 
the H2O ice band, a very volatile ice that cannot exist outside molecular clouds (Figure 11). 



6.1. Need for Laboratory or Theoretical Frequencies for C— H Bonds on the 

Surfaces of Microdiamonds 

Under the assumption that the 3.47 /xm feature is due to C-H bonds, Brooke et al. 
[21] interpreted the correlation with H2O ice as indicating that both C-H bonds and H2O 
ice form in step on molecular cloud dust by hydrogen addition reactions. However, we 
noted that other identifications of the 3.47 ixm feature are also possible. One fertile region 
for astrochemical research is to better establish the formation mechanism and absorption 
spectrum of interstellar "diamond" -like carbon of interstellar size (1 nm; [70]) with surface 
hydrogen atoms. Predictions or measurements of the complete vibrational spectrum, 
including additional modes which could be searched for in astronomical spectra, would be 
particularly useful. It would also be useful to search for alternative identifications for the 
3.47 //m feature with simple molecules, composed of cosmically abundant elements (Fig. 
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1), which would have subhmation temperatures close to that of H2O ice. Again, predictions 
or measurements of other vibrational modes that could be used to test the identification 
would be helpful. 



7. Cg^ in the interstellcir medium 

The fullerene Ceo was first discovered in the laboratory by Kroto et al. [72]. They 
proposed that this novel form of aromatic carbon could potentially play an important role 
in the ISM. Theoretical studies of dust formation in carbon- rich stellar mass-loss [73 - 75] 
suggested that fuUerenes could be formed in such carbon-rich environments and expelled 
into the ISM. Laboratory studies have also shown that one product of the photoerosion of 
hydrogenated amorphous carbon grains is fuUerenes [76]. The existence of fuUerenes such 
as Cgo and C70 in our solar system, for instance in meteoritic samples, is hotly debated [77, 
78]. 

Searches for Ceo in the ISM, through its UV absorption band at 386 nm, have placed 
stringent limits of < 0.01% of the cosmic carbon abundance in Ceo [79, 80]. The dominant 
ionization state of Ceo in the ISM, however, is predicted to be C^q [47, 49, 50]. Foing & 
Ehrenfreund [81, 82] detected diffuse interstellar bands (sec [83] for a review of the diffuse 
interstellar bands) at 958 and 963 nm, which they argued were due to Cgg, based on a 
comparison to laboratory data [84]. 

This inspired Moutou et al. [61] to use ISO data to search for the 7.1 and 7.5 /im 
(1331 and 1406 cm"^) vibrational emission lines of C^q in NGC 7023, a reflection nebula 
with strong lEF emission. Such a search requires high signal-to-noise and high spectral 
resolution (A/AA = 1800), because these features, if present, would appear as weak bumps 
on the blue wing of the strong aromatic 7.7 pLm lEF. From our ISO spectra (Figure 8), we 
place an upper hmit on Cgo in NGC 7023 of < 0.3% of interstellar carbon [61]. 

The central star of NGC 7023, HD 200775, is known to have abnormally weak diffuse 
interstellar bands [85]. If the diffuse interstellar bands observed at 958 and 963 nm are 
due to C6"o, then the general weakness of the diffuse interstellar bands in NGC 7023 might 
naturally lead to a failure to detect Cgg in the mid-infrared toward this source. Astronomers 
can place more robust limits on the abundance of fuUerenes in the ISM by searching for 
the 7.1 and 7.5 /im Cqq lines towards a larger number of lines of sight, particularly those in 
which the diffuse interstellar bands are strong. 
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7.1. Laboratory and Theoretical Work Needed for FuUerenes 

Laboratory astrochemists can also contribute greatly to the search for interstellar 
fuUerenes by further experiments. The assignment of the interstellar bands at 958 and 963 
nm to Cg^ was based on a comparison of gas-phase CgQ in the ISM with CgQ measured in 
inert-gas matrices, with the effects of matrix shifts on the wavelengths being unknown. 
The 7.1 and 7.5 //m Cgg bands were also measured in inert-gas matrices, again introducing 
an unknown amount of matrix shift to the wavelengths. Laboratory measurements of all 
these bands in the gas phase would help immensely to determine whether detections of 
astronomical absorption features near the laboratory wavelengths of C^q prove or disprove 
the existence of Cqq in the ISM. Finally, by symmetry, there should be two longer wavelength 
vibrational modes of Cg"o which were not measured by Fulara et al. [84]. Laboratory 
wavelengths for these two additional modes would open up new avenues for astronomical 
observations of Cqq. 



8. Conclusions 

I have emphasized in this paper a number of puzzling astronomical observations of 
aromatics in the ISM, and related molecules, where understanding of these problems would 
greatly benefit from new experiments in laboratory astrochemistry and progress in chemical 
theory. The study of the ISM is an area where the interaction of astronomy, physics, and 
chemistry has proven particular fruitful in reaching new conclusions about the universe that 
surrounds us, and I strongly encourage astrochemistry groups to tackle these experimental 
and theoretical problems as part of this on-going interaction. 

I would like to thank Tim Brooke for contributions to the text on the 3.47 //m feature. 
I am grateful to Max Bernstein, Jean Chiar and Farid Salama for communicating data in 
advance of publication, and to Jacquie Keane for sharing a preprint before the refereeing 
process was complete. I also thank James Pizagno for speedy assistance with Figure 3. 
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Figure Captions 

Figure 1 — The cosmic abundances of elements from Cox [1], presented as logio of the 
elemental abundance relative to hydrogen versus the atomic number of the element. Only 
elements with abundances greater than 10^^ relative to hydrogen are illustrated. 

Figure 2— The ISOCAM + CVF spectrum of the reflection nebula NGC 7023 from 
Cesarsky at al. [2], plotted as relative flux density versus frequency in cm~^. The spectral 
resolution is v j i^v — 40. Five of the six interstellar emission features (lEFs) can be seen, 
at 790, 890, 1160, 1300, and 1610 cm'^ (6.2, 7.7, 8.6, 11.3, and 12.7 //m). 

Figure 3 — IRAS broad-band photometric observations {jiUed squares) of i?(12/total), 
the ratio of the fiux in the IRAS 12 //m filter to the total infrared flux, plotted against 
Tegf(star), the effective stellar temperature of the star which illuminates each reflection 
nebula, from Sellgren et al. [43]. Error bars are il-o". Upper limits are 3-0". The curves 
arc the predicted values of i?(12/total) from Pizagno et al. (2000, in preparation), derived 
by convolving the UV and visible absorption curves of neutral naphthalene {dashed curve), 
from Salama & AUamandola [45], and ionized naphthalene {solid curve), from F. Salama 
(1999, private communication), with theoretical energy distributions of stars as a function 
of Teff(star). 

Figure 4— Scaled ISOCAM + CVF spectra of three reflection nebulae, vdB 17 (NGC 
1333; dotted line), vdB 59 (NGC 2068; thick solid line), and vdB 133 {thin solid line), 
plotted versus wavelength in //m, from Uchida et al. [44, 46]. The spectral resolution is 
A/AA = 40. The stars which illuminate each reflection nebula respectively have T^s = 
11,000 K for vdB 17, 19,000 K for vdB 59, and 6,800 + 12,000 K for vdB 133 (which is 
illuminated by a binary star). 

Figure 5 — Ratio of the 5.5 - 9.75 /im fiux to the 10.25 - 14.0 /im flux versus the UV 
radiation field, Gq, from Uchida et al. [46]. All data are derived from ISOCAM + CVF 
spectra. Different nebular positions are illustrated for vdB 17 {filled diamonds; [46]), and 
vdB 59 {filled circles; [46]). Results at a single nebular position are shown for vdB 133 
{open circle; [46]), NGC 7023 {cross; from Cesarsky et al. [2]), and p Oph {open square; 
from Boulanger et al. [15]). Error bars are il-u. 



Figure 6 — ISOCAM + CVF spectra of two reflection nebulae, plotted versus 
wavelength in /xm, from Uchida et al. [46]. The right panel is vdB 59 (NGC 2068), while 
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the left panel is vdB 17 (NGC 1333). Spectra at different nebular positions are labeled 
by their values of the UV intensity, Gq. The spectral resolution is A/AA = 40. Note the 
broadening of the 7.7 /im lEF at the lowest Gq values in vdB 17, which is unexplained by 
current models. 

Figure 7 — FWHM in fim of the 6.2, 8.6, and 7.7 fim lEFs as a function of the UV 
radiation field. Go, within vdB 17, from Uchida et al. [46]. Data are derived from ISOCAM 
+ CVF spectra. Error bars are ±1-(t. No correction has been made for the instrinsic 
spectral resolution of the instrument (A/AA = 40). 

Figure 8 — Spectra at 7.0 - 8.7 //m of NGC 7023, plotted versus wavelength in /xm, 
observed with ISO's Short Wavelength Spectrometer (SWS) at A/AA = 1800 by Moutou 
et al. [61]. Solid vertical fines mark the laboratory positions of Cgg bands [84]. The data 
labeled "(up - down)/2" provide an estimate of the noise of the spectrum. The emission line 
in NGC 7023 at 8.02 /xm is tfie 0-0 S(4) line of molecular fiydrogen. The SWS responsivity 
calibration is shown at the bottom of the plot, offset for clarity. 

Figure 9 — (top) Optical depth, plotted versus wavelength in fim, of a possible 
aromatic absorption feature toward the protostar Mon R2/IRS-3, embedded in a molecular 
cloud, from ground-based spectra (A/AA = 1000) obtained by Sellgren et al. [19] at 
the United Kingdom Infrared Telescope Facility, (middle) Optical depth, plotted versus 
wavelength in /im, of a possible aromatic absorption feature toward the enigmatic Galactic 
Center source GCS3, which is obscured primarily by the diffuse interstellar medium, from 
Short Wavelength Spectrometer spectra (A/AA = 1000) from ISO observed by Chiar et al. 
[66]. [bottom) Normalized intensity of the 3.3 ^ra lEF toward the refiection nebula NGC 
7023 (Sellgren et al. 2000, in preparation), from tfie Sfiort Wavelengtfi Spectrometer on ISO 
(A/AA = 1000). Tfie emission line at 3.234 /xm is tfie 1-0 0(5) line of molecular fiydrogen. 

Figure 10 — Tfie optical deptfi of tfie 3.47 fira absorption band in molecular clouds, 
possibly due to "diamond" -like carbon, plotted versus tfie optical deptfi of 9.7 /xm silicate 
absorption, a refractory grain component. Data are from Brooke et al. [22] {filled circles), 
Sellgren et al. and Brooke et al. [18, 21] {open circles), and Cfiiar et al. [82] {crosses). 

Figure 11 — Tfie optical deptfi of tfie 3.47 iim absorption band in molecular clouds, 
possibly due to "diamond" -like carbon, plotted versus tfie optical deptfi of 3.1 /xm H2O ice 
absorption, a volatile grain component. Data are from Brooke et al. [22] {filled circles), 
Sellgren et al. and Brooke et al. [18, 21] {open circles), and Cfiiar et al. [82] {crosses). 
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